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Abstract 

We consider optimal investment problems for a diffusion market model with non- 
observable random drifts that evolve as an Ito's process. Admissible strategies do not use 
direct observations of the market parameters, but rather use historical stock prices. For a 
non-linear problem with a general performance criterion, the optimal portfolio strategy is 
expressed via the solution of a scalar minimization problem and a linear parabolic equation 



> 

with coefficients generated by the Kalman filter. 
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The paper investigates an optimal investment problem for a market which consists of a locally 
risk free asset, bond or bank account with interest rate r(t), and a finite number, re, of risky 
stocks. We assume that the vector of stock prices S(t) evolves according to an ltd stochastic 
differential equation dSi(t) = Si(t)[a,i(t) dt + J2j a ij(t) dwj(t)], i = l,...,re, with a vector 
of appreciation rates a(t) and a volatility matrix cr(t). The problem goes back to Merton 
(1969), who found strategies which solve the optimization problem in which EU(X(T)) is to be 
maximized, where X(T) represents the wealth at the final time T and where U(-) is a utility 
function. If the market parameters are observed, then the optimal strategies (i.e. current vector 
of stock holdings) are functions of the current vector (r(t),a(t),a(t),S(t),X(t)) (see, e.g., the 
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survey in Hakansson (1997) and Karatzas and Shreve (1998)). But in practice, a(t) and a(t)) 
have to be estimated from historical stock prices or some other observation process. There are 
many papers devoted to estimation of (a(-), cr(-)), mainly based on modifications of Kalman- 
Bucy filtering or the maximum likelihood principle (see e.g. Lo (1988), Chen and Scott (1993), 
Pearson and Sun (1994)). Unfortunately, the process a(-) is usually hard to estimate in real-time 
markets, because the drift term, a(-), is usually overshadowed by the diffusion term, cr(-). On 
the other hand, a(t) can, in principle, be found from stock prices. Thus, there remains the 
problem of optimal investment with unobservable a(-). 

In fact, the problem is one of linear filtering. If Ri(t) is the return on the ith stock, then 
dR(t) = a(t)dt + a(t)dw(t), so the estimation of a(t) given {R(t),t < t} (or {S(t),t < t}) is 
a linear filtering problem. If a(-) is conditionally Gaussian, then the Kalman filter provides the 
estimate which minimizes the error in the mean square sense. 

A popular tool in optimal control and filtering theory is the separation theorem. This the- 
orem has an analog in portfolio theory: it is the so-called "certainty equivalence principle": 
agents who know the solution of the optimal investment problem for the case of directly ob- 
servable a(t) can solve the problem with unobservable a(t) by substituting E{a(t)|5(r), r < t} 
(see, e.g., Gennotte (1986)). Unfortunately, this principle does not hold in the general case of 
non-log utilities (see Kuwana (1995)). Note that this principle is unrelated to the much more 
recent notion of "certainty equivalent value" to be found in the work of Frittelli (2000). 

Williams (1977), Detemple (1986), Dothan and Feldman (1986), Gennotte (1986), Brennan 
(1998) solved the investment problem using the Kalman-Bucy filter and dynamic programming. 
By this method, the optimal strategy can be calculated via solution of the Bellman parabolic 
equation; this equation is non-linear. 

Karatzas (1997), Karatzas and Zhao (1998), Dokuchaev and Zhou (2000), Dokuchaev and 
Teo (2000) have obtained optimal portfolio strategies in general non-Gaussian setting, but only 
for case of time independent coefficients. 

An approach based on Malliavin calculus gives a possibility to consider more general setting. 
Lakner (1995), (1998) assumes that S(-) and w(-) have equal dimension (as we do), and that r(-) 
and cj(-) are deterministic. This again guarantees that the filtration of S(-) is Brownian. Results 
from filtering theory give a representation of the optimal portfolio, which is explicit in terms 
of a conditional expectation of a Malliavin derivative when the cij(-) are Ornstein-Uhlenbeck 
processes independent of w(-). Karatzas and Xue (1990) assume that there are more Brownian 
motions than stocks. They assume that r(-) and cr(-) are adapted to the observable S(-). After 
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projecting onto an n-dimensional Brownian motion which generates the same filtration as S(-), 
they obtain a reduced, completely observable model; existence of an optimal portfolio follows, 
but the optimal strategy is, as usual, defined only implicitly. 

We also consider the optimal investment problem with random and unobservable a(-). Fol- 
lowing Lakner (1998) and Rishel (1999), we assume that a(t) is a Gaussian process modelled by 
a system of linear Ito's equations. However, we consider a more general case when (a(-),r(-)) 
may depend on the realized returns (i.e., b(-) ^ in equation (|2.4p below, and r(-) is correlated 
with £(•)). We express the optimal strategy via solution of a Cauchy problem (|4.3|) . (|4.8|) for 
a linear parabolic equation in (n + l)-dimensional vector space. Thus, we propose a simpler 
method than dynamic programming: the nonlinear parabolic Bellman equation is replaced for a 
linear parabolic equation. Note that the solution in Lakner (1998) expresses the optimal strategy 
via a conditional expectation of a random claim that depends on w(-); the solution presented 
below is also based on the martingale method but it is more constructive provided we can solve 
the Cauchy problem (|4,3p . (|4,8p . Using the technique of backward stochastic partial differential 
equations, we prove existence and uniqueness of the solution for this Cauchy problem. Further- 
more, the most restrictive condition in Lakner (1998) was that the initial covariance of a(0) is 
small enough (the condition (13. 5p below). We replace it by another condition (|4.9p that depends 
on U: it is less restrictive than (]3.5f) for some C/'s and more restrictive for others C/'s. For some 
problems, our condition (|4.9p is automatically satisfied. In addition, we allow correlated a(-) 
and w(-). 



Consider a diffusion model of a market consisting of a locally risk free bank account or bond 
with price B(t), t > 0, and n risky stocks with prices Si(t), t > 0, i = 1,2, ... ,n, where n < +co 
is given. The prices of the stocks evolve according to the following equations: 



where Wi(t) are standard independent Wiener processes, aj(i) are appreciation rates, and aij(t) 
are volatility coefficients. The initial price Si(Q) > is a given non-random constant. The price 
of the bond evolves according to the following equation 



2 The Model and Definitions 




t > o 



(2.1) 




(2.2) 
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where B(0) is a given constant which we take to be 1 without loss of generality, and r(t) is the 
random process of the risk-free interest rate. 

We are given a probability space (fi, T , P), where Q is a set of elementary events, T is a 
complete cr-algebra of events, and P is a probability measure. 

We introduce the vector processes ( T denoted transpose) 

w(t) = {w l {t),... J w n {t)) T , S(t) = (Si (£),..., <S„(t)) T , a(t) = (ai(t),...,a n (t)) T , 

and the matrix process a(t) = {<Tij(t)}™j =1 . 

Let 1 = (1, ... , 1) T G R n , and a(t) = a(t) - r(t)l. 

We define the return to time t by dRi(t) = dSi(t)/ Si(t), Ri(0) = 0, and introduce the vector 
of returns R(t) = (R\(t), . . . , R n (t)) T and of excess returns Ri(t) = Ri(t) — J Q * r(s) ds. 

Let {J r t' r }o<t<T be the filtration generated by the process (r(t),S(t)) completed with the 
null sets of T. 

Set S{t) = exp (- f* r(s)da) S(t). 

We denote by |x| the Euclidean norm of a vector x € R fc . For an Euclidean space E, we 
denote by B([0,T}; E) the set of bounded measurable functions f(t) : [0,T] — > £7. We denote 
by /„ the identity matrix in R nxn . As usual, we say that A < B for symmetric matrices if the 
matrix B — A is definitely positive. We denote <fi~ = max(0, — (f>), and we denote by the 
indicator function. 

The model for r, a, and a 

To describe the distribution of ait), we shall use the model introduced in Lakner (1998, p. 84), 
generalized for our case of random r, non-constant coefficients for the equation for a, and 
correlated r, a, and w. We assume that we are given measurable deterministic processes a(t), 
0(t), b(t) and 5(t) such that 

da(t) = a{t)[S(t) -a{t)]dt + b(t)dR(t) + (3(t)dWit), (2.3) 

where a(t) G R nxn , 0{t) G R nxn , b{t) G R nxn , 5(t) G R n , and where is an n-dimensional 
Wiener process in (Q, P). We assume that a{t), Pit), bit), and 5(t) are continuous in t and 
such that the matrix (3(t) is invertible and |/3(t) _1 | < c, where c > is a constant. Further, 
we assume that a(0) follows an n-dimensional normal distribution with mean vector niQ and 
covariance matrix 70. The vector ijiq and the matrix 70 are assumed to be known. We note 
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that this setting covers the case when a is an n-dimensional Ornstein-Uhlenbeck process with 
mean-reverting drift. 

Clearly, equation (|2.3|) can be rewritten as 

da(t) = (a(t)8(t) + [b(t) - a(t)]a(t)) dt + b(t)a(t)dw(t) + (3{t)dW(t). (2.4) 

In addition, it can be seen that Ri(t) evolve as 

n 

dRi(t) = Oi(t)dt + ^2<T ij (t)dwj(t), t > 0. (2.5) 
j=i 

We assume that the process a(t) is continuous in t, non-random and such that a(t)a(t) T > 
c a I n , where c CT > is a constant. 

Further, we assume that r(-) = <p r (R(-),Q), where G is a random element in a metric space 
X r , and where <p r : C([0,T];R n ) x X T — > B([0, T];R) is a measurable function, and does not 
depend on (w(-),W(-),a(0)). In addition, we assume that the process r(t) is adapted to the 
filtration generated by (R(t),@). Note that closed system (|2.4j) - (j2.5|) for the pair (a(t),R(t)) 
does not include r(-), and (a(-),R(-)) does not depend on 0. Therefore, the market model is 
well defined. The assumptions for measurability of r don't look very natural. However, they 
cover generic models when r is independent on R or non-random, and we can still consider some 
models with correlated r and R. 

Under these assumptions, the solution of (|2.ip is well defined, but the market is incomplete. 

Let (p m {t, s), m = 0, 1, be the solution of the matrix equation 

%(i, s) = [m ■ b(t) - a(t)]4> m {t, s), 

<t>m{s,s) = I n . 

Let 

Km(t)= [ 0m(t,s)b(s)a(s)a(s) T b( S ) T m (t,s) T ds, m = 0, 1. (2.6) 
J o 

We have that 

a(t) = Mt, 0)5(0) + I 0i (i, s)[a{s)5{s)ds + b(s)a(s)dw(s) + /3(s)dW(s)}. 
Jo 

It follows that K\(t) is the covariance matrix for a(t) calculated with (3(t) = and a(0) = 
0. By the linearity of (|2.4p . it follows that K\(t) is the conditional covariance for a(t) given 
(^(■)| [0 ,t],a(0))or (W(-)\ [0:T] ,a(0)). 

Note that K m (t) can be found as solutions of linear equations that one can easy derive from 
dZ3D and flU]) (see, e.g., Arnold (1973), Chapter 8). 
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We assume that b is "small". More precisely, we assume that there exists e > such that 
TK m {t) +el n < a(t)a(t) T Vie[0,T], m = 0, 1. (2.7) 

The risk neutral probability measure 

Set Q(t) = ( ( j(t)fT(i) T )- 1 , and set 

Z = exp (J [a(t)' 1 a(t)] T dw(t) + 1 jf a(t) T Q(t)a(t)d?j . (2.8) 

Proposition 2.1 

Ejexp iy a(t) T Q(t)a(t)dt W(-),a(0)} < +oo a.s. (2.9) 

By this proposition, the Novikov's condition is satisfied conditionally, E{Z~ X | W(-), a(0)} = 
1, then EZ^ 1 = 1. 

Define the (equivalent) probability measure P* by dP*/(iP = Z~ x . Let E* be the corre- 
sponding expectation. 

The wealth and strategies 

Let Xq > be the initial wealth at time t = 0, and let X(t) be the wealth at time t > 0, 
X(0) = Xq. We assume that 

n 

X(t)=7ro(t)+J>i(t), (2-10) 

i=l 

where the pair (7To(£),7r(t)) describes the portfolio at time t. The process ttq{€) is the investment 
in the bond, Hi(t) is the investment in the ith stock, ir(t) = (vri(t), . . . ,ir n (t)) T , t > 0. 

Definition 2.1 The process X(t) = exp ^— J * r(s)dsj X(t) is called the normalized (or dis- 
counted) wealth. 

Let S(t) = diag (S\(t), . . . , S n (t)) and S(t) = diag (Si(t), . . . , S n (t)) be diagonal matrices 
with the corresponding diagonal elements. The portfolio is said to be self-financing, if 

dX{t) = 7r(t) T S(t) -1 dS(t) + 7r (t)r(t)dt = Tr{t) T dR(t) + TT (t)r(t)dt. (2.11) 

It follows from (|2.1U|) that for such portfolios 

dX(t) = r(t)X(t) dt + Tr(t) T (a(t) dt + a(t) dw(t)) , ^ 
dX(t) = B(t)~ l 7r(t) T dR(t), 
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so 7r alone suffices to specify the portfolio; the process ttq is uniquely defined by ir via 
(|2.10p . (|2.12p ; 7r it is called a self-financing strategy. 

Definition 2.2 Let £ 6e the class of all J- ~ t ' T -predictable processes ir(-) such that 

• lo (l 71 "^) 7 "^)! 2 + k(t) T cr(t)| 2 ) dt < oo a.s. 

• there exists a constant q n such that P (x(t) — Xq > q n ,\/t E [0,T]^ = 1. 

A process 7r(-) € S is said to be an admissible strategy with corresponding wealth X(-). 

For an admissible strategy tt(-), X(t,Tr(-)) denotes the corresponding total wealth, and 
X(t, vr(-)) the corresponding normalized total wealth. It follows that X(t,ir(-)) is a P*- 
supermartingale with E*X(t, vr(-)) < X and E*|X(i, vr(-))| < |X | + 2\q n \. 

Note that by definition, admissible strategies from S use observations of r(t) and only. 
For these strategies, the processes X(t) and X(t) are F t ' r -adapted. 

The following definition is standard. 

Definition 2.3 Lei ^ be a given random variable. An admissible strategy ir(-) is said to replicate 
the claim £ if X(T,n(-)) = £ a.s. 



3 Problem statement and preliminary results 

Let T > 0, let D C R be convex and bounded below, and let Xq G D be given. Let [/(•) : D — > 
R U {— oo} be such that U(Xq) > — oo. 

We may state our general problem as follows: Find an admissible self-financing strategy ir(-) 
which solves the following optimization problem: 

Maximize EU(X(T, tt(-))) over vr(-) G £ (3.1) 

( X(0, ?r(-)) = X , 
subject to < „ ^ (3-2) 

[ X(T,tt(-)) £ D a.s. 

The condition X(T, 7r(-)) € -D may represent a requirement for a minimal normalized termi- 
nal wealth if D = [k, +oo), fc > 0. This condition may represent also a requirement for the 
normalized terminal wealth in goal achieving problems if D = [ko, k±], ko < k\. 
We assume that U, Xq and D satisfy the following two conditions. 
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Condition 3.1 There exists a measurable set A C [0, oo), and a measurable function F(-,-) : 
(0, oo) x A — > D such that for each z > 0, x = F(z, A) is a solution of the optimization problem 

Maximize zU{x) — Ax over x G D. (3-3) 

Note that the usual concavity hypotheses imply this condition, but more general utility functions 
are also covered. For example, this condition is satisfied for the goal achieving problem when 
U(x) is a step function (see e.g. Karatzas (1997), Dokuchaev and Zhou (2000)). 

Let Z = ~E{Z\Fj,' r }. Since (R(-),a(-)) does not depend on G, we have that Z does not 
depend on 0, and Z = E{Z\R(-)}. Let F(-) be as in Condition 13.11 

Condition 3.2 There exists A G A such that E*\F(Z, A)| < +oo and E*F(Z, A) = X . 

We solve our problem in two steps using the martingale approach. First we show that 
EU(F(Z, A)) is an upper bound for the expected utility of normalized terminal wealth for 
7r(-) € S. Then we find a portfolio 7r(-) which replicates the claim B(T)F(Z, A). This establishes 
the optimality of 7?(-). 

The optimal claim 

The following theorem is a reformulation of Theorem 2.5 from Lakner (1998) under slightly more 
general conditions that allow discontinuous functions F and U such as step functions. 

Theorem 3.1 (Dokuchaev and Haussmann (2000)). With A as in Condition \3.2\ let £ = 
F(Z,X). Then 

(i) E [/-(£) < oo, £ G 5 a.s.; 

ft) Etf(0 > EC/(X(T,7r(-))), Vvr(-) G S; 

^Mij T/ie claim B(T)^ is attainable in E, and t/iere exists a replicating strategy in E. T/iis 
strategy is optimal for problem li 3. 1\) 

This theorem uses duality approach for constrained optimization that goes back to Lagrange, 
and A is the corresponding Lagrange multiplier. 

Remark 3.1 Theorem 2.5 from Lakner (1998) was stated under some additional assumptions 
that can be formulated in our notations as 

(i) b(t) = 0, r is non-random, r, a, a, f3, 5 are constant, and D = (0, +oo); 
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(ii) U is strictly concave and continuously differentiable on (0, +00), and lim x ^ +00 U'(x) = 0; 

(iii) there exists a function J(-) : D — > R such that J(X/x) = F(x, A); 

(iv) E„J(X/Z) < +00 for any A > 0. 

Solution via conditional expectation 

Let 

o(t) ± E{a(t) | p }- 
Set 5(0 = a(t) — b(t) and tuq = Ea(0). 

Let j(t) £ R, nXfl j-jg ^he unique solution (in the class of symmetric nonnegative definite 
matrices) of the deterministic Riccati's equation 

g(t) = -[&(t)ff(t) T + 7 (t)]Q(*)[&(*M*) T +7W] T " 5(t) 7 (t) -7(t)5(t) T + /3(W) T , 
7(0) = 70. 

(3.4) 

Here 70 = E[5(0) - m ] [5(0) - m ] T . In fact, j{t) = e{[5(0 - a(t)][a(t) - a(t)] T |J c f' r }. 
Let A(t) = —5(0 — j(t)Q(t), and let <fi(t) be the solution of the matrix equation 

f t (t) = Aimtl 

x 0(0) = /n, 

where I n is the unit matrix in R nxn . 

The following theorem is a reformulation of Theorem 4.3 from Lakner (1998). It gives the 
solution of the investment problem via conditional expectation of future values of some processes 
with known evolution. 

Theorem 3.2 (Lakner (1998)). Let conditions (i)-(iv) in Remark \3. 1\ holds, let U(x) be twice 
differentiable on (0, +00) , and let 

tr 7 o + T||/3|| 2 <^ 1 , K 1= 1 , K = max ||e- a T, (3.5) 

360T||cr 1 |p-Ko *g[o,t] 

where \\ ■ \\ denotes the Frobenius matrix norm, i.e., ||c _1 || 2 = tr [cr~ • L cr -lT ]. Further, let 

J{x) < K(l + x~ 5 ), - J'(x) < K(l - x' 2 ) (3.6) 
for some K > 0. Then the optimal strategy is 
7r(t) T = £T(t)Z(t)E j J'(\Z)Z~ 2 

where H(t) = \e r ^~ T ^Q and w{t) = w(t) — /* a~ 1 a(s)ds. 



--(ON/,) 1 ]- 1 Vr 1 ^*) -2(0 



■jrS,r 

t (1 



We propose below another solution such that the optimal strategy is presented via solution 
of a linear deterministic parabolic equation. We replace conditions (|3,5p by condition (|4,9p that 
can be less restrictive and is always satisfied if D is bounded. In addition, we dropped condition 
()3.6p and the condition that (r, a) and w are independent: we allow b(-) ^ and r = 4> r (R(-), G). 



4 Main results: Solution via linear parabolic equation 

Let y(t) = (yi(i), . . . ,y n+ i(t)) = (a(t),y n+ i(t)) be a process in R n+1 , where 

a(t) = E&{t)\jf r }, 
Vn+i{t) = ~ f a(s) T Q(s)a(s)ds + [ a(s) T Q{s) dR(s). 



Jo 



Let functions /(•) : R n+1 x [0,T] R n+1 and g(-) : R™ +1 x [0,T] R(^+i)x« be such that 

ix T Q(t)£ J ' ^ X ' ^ x T Q(t) 

Here and 7(4) are matrices defined above, 7(i) is the solution of (|3.4p . and 

x = (xi, . . . ,a? n+ i) 1 = I ), x = (xi, . . . ,x n ) T . 

X n +l 

By Theorem 10.3 from Liptser and Shiryaev (2000), p. 396, the equation for a(t) is 



da(t) = [A(t)a(t) - b(t)a(t) T Q(t)a(t) + a(t)5(t)]dt + [b(t)a(t) T + j(t)]Q(t)dR(t), 
a(0) = niQ. 

By (I4.ip - (|4.6I) . it follows that y(-) is the solution of the Ito's equation 

dy(t) = f(y(t),t)dt + g(y(t),t) dR(t), 
y(o) = yo, 

with 



(4.1) 



(4.2) 



Vo= I ^° j €R n+1 , m = E5(0). 

The function f(y,t) here does not satisfy Lipshitz condition with respect to y £ R n+1 . However, 
the solution of this equation is uniquely defined. (It is shown in the proof of Lemma 14.11 below 
that the solution of (|4.2p can be presented as a part of the unique solution of some Ito's equation 
with coefficients that are affine with respect to the state variable). 
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Lemma 4.1 Let a function $(•) : R n+1 — > R be such that 
(i) VMV(T)) = X ; 

(ii) is continuously twice differentiable; 

(Hi) E^(y(T)) 2 < +oc. 

Then there exists a unique classical solution V : R n+1 x [0, T] — > R of the boundary value 
problem 

dV dV 1 f d 2 V 1 

-Qjrfa t) + — (x, t)f(x, t) + -tr I -q^{x, t) g(x, t)a(t)a(t) T g(x, t) T J = 0, (4.3) 

V(x,T) = $(x). (4.4) 

Further, the processes X(t, 7r(-)) = V(y(t),i) and vr(t) T = B(t)^(y(t),t)g{y(t),t) , are uniquely 
defined as elements of the spaces C([0, T], L2(f2, .F, P*)) and ^([0, T], ^(n, P, P*)) respectively, 
and there exists a constant C > suc/t i/iai 

sup E,|X(t,vr(-))| 2 + E* / B(t)~ 2 \Tr(t)\ 2 dt < CE*|$(y(T))| 2 (4.5) 
te[o,T] Jo 

for all these <£. Furthermore, the strategy ir(t) = (wi(t), . . . ,ir n (t)) belongs to £ and replicates 
the claim B(T)&(y(T)) given the initial wealth X$ with the normalized wealth X{t) = V(y(t),t). 

Note that estimate (|4.5p reminds the Krylov-Ficera estimate (see Theorem 5.3.3 from Ro- 
zovskii (1980)) or its modification from Dokuchaev (1995)). 
Further, we have that 

dZ(t) = a{t)Z{t)dR{t). (4.6) 

This formula (j4.6l) was derived in Theorem 3.1 from Lakner (1998) for the case when a is 
constant and 6 = 0. The proof for a non-constant o~(t) and b ^ can be found in Dokuchaev 
and Haussmann (2000) and in Chapter 9 from Dokuchaev (2002). It follows that 

y n+ i(t) = ]nZ(t). (4.7) 

Introduce the function e(-) : R n+1 — ► R such that e(y) = exp[y n+ i] for y = (y\, . . . ,y n+ i) T . 
Note that Z = e(y(T)). 

Let V = V(x, t, A) : R n+1 x [0, T] x A — > R be the solution of partial differential equation 
(14. 3D with the condition 

V(x,T,\)=F(e(x),\). (4.8) 
The following result now is immediate. 
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Theorem 4.1 Let X be such as in Condition \3.2l Assume that the function F(-, A) : R — > R is 
suc/i i/iai is suc/i i/iai conditions (i)-(ii) of Lemma \4-1\ are satisfied with <£(x) = F(e(x),X), and 

E,F(J,A) 2 < +00. (4.9) 

Then there exists an unique classical solution V of problem f or A = A, and there 

exists an admissible self-financing strategy ir(-) G £ which replicates the claim B(T)F(Z, A). 
TTiis strategy is an optimal solution of problem f 3. 1\) and 

n(t) T = B(t)^-(y(t),t, X)b(y(t),t), X(t, vr(-)) = V(y(t),t, A), (4.10) 

Note that it is possible that condition (|4.9p is not satisfied but the optimal claim F(Z, A) is 
still replicable in the class of strategies S. For example, let U(x) = logx, Xq = 1, and (0, +00) C 
D, then A = (0,oo), F(z, A) = z/X, X = 1, and the strategy is vr(t) T = B(t)a(t) T Z(t)Q(t) is 
replicating (and optimal) even in the case when (|4,9p is not satisfied. 



5 Special cases 

Note that conditions (13. 5h were imposed in Lakner (1998) with the only purpose to ensure that 

E*Z 5 < +00, E*Z- 4 < +00. (5.1) 

Our condition (|4.9p for examples (i)-(iii) listed below is satisfied if < +00 for some 

/i G R. For example (i), condition (|4.9p is less restrictive than (|5.ip if I < 5/2 and more 
restrictive if I > 5/2. For example (ii), condition (|4.9p is less restrictive than (|5.ip if Z < 2 and 
more restrictive if / > 2. For example (iii), condition (|4.9p is always less restrictive than (|5.ip . 
These examples are from Dokuchaev and Haussmann (2001): 

(i) Power utility. Assume D = [0, +00), Xq > 0, U(x) = d~ l x d , where either d G (0,1) or 
d < 0. Then A = (0, 00), F(z, A) = (z/X) 1 , and A = X Q 1/l (E^Z 1 ) 1 / 1 , where Z = 1/(1 - d). 

(ii) Assume D = [0, +00), U(x) = —x d + x, where d = 1 + l/l, and Z > is an integer, 
X > dT l . Then A = [0, 00), F(z, A) = (1 + X/z) l d~ l , X is a root of a polynomial of degree Z. 

(iii) Mean-variance utility. Assume D = R, ZJ(x) = — /err 2 + cx, where k G R and c > 0, 
X > 0, then F(z,X) = (c - X/z)/(2k). 

We present below some sufficient conditions that ensure E*Z^ < +00 and, therefore, can be 
useful for verifying (16. 2p . 

Let K(t) be the covariance for a(t) under the probability measure P*, and let K(t) be the 
covariance for a(t) under P*. 
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Lemma 5.1 If p G [0, 1], then E*iJ M < +oo. Let p, < or p, > 1. T/ien E*iJ M < +oo f/ t/iere 
exist e > and p > 1 suc/i t/iai a£ /east one of the following conditions holds: 

(i) K{p)K(t) < a(t)a(t) T -el n fort G [0,7], where n{p) = qT{p 2 p-p) > with q = p(p-l) _1 . 

(iij K(p)K(t) < a(t)a(t) T - el n for t G [0,T]. 

It follows from Proposition 17.21 below that K(t) and if(t) are the covariances of the processes 
defined by (|2.4p and (|4,ip respectively with R(-) replaced by i?*(-). Thus, these covariances can 
be found as solutions of linear deterministic equations that one can easy derive from (|2.4p and 
flU} (see, e.g., Arnold (1973), Chapter 8). 



6 Case of discontinuous F 

To proceed further, we shall need a special weighted L2-space with a weight defined via some 
parabolic equation. First, we introduce the operator 



n+1 q ^ n+1 q2 

M(t)p t - £ _ (p(x)/i(x,t)) + - ^ 5^ ( P (x)g i3 (x,t)) , 

1,7 = 1 * J 



i=l i,j 

where g = go~a T g T . 

Let ^ G L 2 (R n+1 ) n C 2 (R n+1 ), i = 1,2, be given such that pi{x) > for all x G R n+1 and 
Jr»+i Pi(x)dx = 1. 

We consider the following parabolic equation 

fE(s,i)=A4(t)p(:r,i)+pi(z), tG[0,T], 
p(x,0) = pofc)- 

This boundary value problem has the unique classical solution p(x, t) that is continuous in 
R n+1 x [0,T]. Let 

p(x) = min p(x, t). 

t€[0,T] 

We have that 

p(;t) = G(t,Q)p + [ G(t,s) Pl ds, 
Jo 

where G(t, s) is the semigroup operator generated by (16. ip (with pi = 0) and such that 
G(s,s)pi = pi. We have that (G(t, s)pi)(x) > for i G [s,s + e) for some e = e(x,s) > 0. 
Hence p(x,t) > for all and p(x) > for all x G R n+ . We shall use this p as a weight 
function. 
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We have that p G L 2 (R n+1 ) n Li(R n+1 ), since \p(x)\ < \p (x)\. 
We introduce the weighted space L2,p(R n+1 ) with the norm 

1/2 



We introduce the space 34 of functions u = {ui(x, t)}f =1 : R n+1 x [0, T] — > R fc with the norm 



\=1 



k p J> 








1/2 



Further, we introduce the space W 1 of functions u = u(x, t) : R n+1 x [0, T] -> R with the norm 



<9x ^ 



Finally, we introduce the space Wq consisting of all functions u(-) S W 1 such that u(-) S 
C([0,r];L 2 , P (R n+1 )) with the norm 



w i= sup |H-, *)|| L2ip(Rn +i) + IMIw 1 - 



A 

tG[0,T] 



The above space is a Banach space, since the weighted space L2 !(3 (R n+1 ) is a Hilbert space. 

In fact, the spaces 34 > W , and W^, are the completions in the corresponding norms of the 
set of smooth functions u : R n+1 x [0, T] — > R fc or u : R™ +1 x [0, T] — > R respectively that have 
finite support. 

Theorem 6.1 Let p be the solution of jl6.1\) . and let be the corrsponding space defined via 
the weight p(x) = min^^T] p(x, t). Let $(•) : R n+1 — > R be a measurable function such that 

p{x,T)<$>(x) 2 dx < +oo. (6.2) 



/Rn + 1 

Then boundary value problem ^4-3\ )- [4-4\ ) admits a unique solution V £ W^. Moreover, there 
exists a constant C > independent on <!>(•) and such that 

\\V\\ 2 wl <c[ p(x,T)<S>(x) 2 dx. (6.3) 



'a 



Rn+1 



Note that condition (16.21) allows discontinuous <£. 



Remark 6.1 The definition of Wq ensures that problem (|4.3p - (|4.4p can be stated in Wq. The 
functions V and (dV/dx)g are measurable and Z^p-integrable. The equality in (14. 4p is the 
equality for elements of the space L2 )P (R n+1 ), it is meaningful since V(-,t) is continuous in t in 
L 2 ,p(R n+1 ). The equality in (14.31) is the equality for elements of the dual space W 1 *, since all 
components of ^r(x,t) g(x,t)a(t)a(t) T g(x,t) T belong to W 1 *. 
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It follows from the proof of Theorem 16.11 below that \\p(-, T) ||l 1 (r™+i) = 2. Hence (|6.2p 
is satisfied for any bounded <I>. In addition, it can be shown that \\p(-, T)\\ L2 (^ Rn +i^ < 
CSi=i2 ll/ 9 «(')llL2(R ri+1 )' wnere C > is a constant that does not depend on p^. Therefore, 
(p^D is satisfied for any 3> e L 4 (R n+1 ). 

Theorem 16.11 gives the possibility to present the optimal investment strategy via solution of 
(|4.3p ~- (|4.4p for the case of discontinuous F. An example is the goal-achieving problem, when 
D = [0,06), X 6 (0,a), and U(x) = if < x < a, U{x) = 1 if x > a. Then A = (0,oo), 
F(z, A) = a if < A < z/a, F(z, A) = if A > z/a, and §672$) holds for = F(e(x), A) (VA). 



7 Appendix: Proofs 

Proof of Proposition \2.1\ . By Jensen' inequality, it follows that 

E{expi/ T H(t) T Q(t)5(t)dt W{-), a(0)} = Ejexp i ^a(t) T Q(t)a(t)dt W(-),a(0)} 
< i J Q T Ejexp ^a(t) T Q(t)a(t)dt W(-),a{Q)\. 

We have for definitely positive matrices that if A > i? > then > y4 _1 . By condition (|2.7p 
with m = 1, it follows that 

ifo(t)- 1 > T[a(t)a(t) T - sin]- 1 = TQ(t)[I n - eQ(i)]- 1 

= TQ(i) [/„ + ESWfi)}"] > rO(t) + TeQ(t) 2 > TQ(t) + M, 

where M = M(e) > is a definitely positive constant matrix. Clearly, we can take e > small 
enough to ensure convergency of the series in (|7.ip . 

To complete the proof, we shall use the following fact. Let £ be a Gaussian n-dimensional 
vector, Kg = E(£ — E£)(£ — E£) T > 0. It is known that the probability density function for £ is 
Cexp[-i(x-E£) T ^ 1 (o;-E£)], where C > is a constant. It follows that E exp(i£ T P£) < +oo 
for any matrix P £ j^ nxn such that < P < K^ 1 . Then the proof follows from (|7.ip , □ 

We introduce the process 

R t (t) = [ a(s)dw(s). 



o 



Let n-dimensional vector random process a*(t) be defined as the solution of 

da,(t) = (a(t)6(t) - a(t)a*(t))dt + b(t)dR,(t) + (3(t)dW{t), 2,(0) =5(0). 

Set T T 

Z* = exp (J [(7(t) _1 o,(t)] T du;(t) - ^ jf S*(t) T Q(i)S*(t)d^ . (7.2) 
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Proposition 7.1 There exists a measurable function ip : C([0,T];R n )xB([0,T];R n ) — ► R such 
that Z* = ij)(R*(-), a*(-)) and Z = ip{R{-),a{-)). 

Proof. Clearly, ip is denned by 

logZ* = j£ a*(t) T Q(t) - ia*(t)d^ . (7.3) 

□ 

Let r*(-) = 4> r (R lf {-), 0) and -B*(i) = l?(0)exp ^/ * ^(s)(is^ (0 r is denned in Section [2]). Let 

= E{Z*\R*(-),r *(■)}. (7.4) 

Let T = C([0, T];R n ) x R n . Clearly, there exists a measurable mapping „4 : [0, T] x 
C([0,T];R n ) x T -> C([0,T];R n ) such that S»(t) = .4(t, #*(•), W(-), 2(0)) and o(t) = 
^(t,i?(-),^(-),5(0)). 

We have that Z* = E{Z*\R*(-)} = and 

Z* = E{iP([R*(-),a*(-)]\R*(-)} = Ety[iL(-),A(-,R,(-),W(-)Mfy)]\M-)}' 
By Proposition I7-H it follows that 

Z = E{ip[R(-),a(-)]\R(.)} = E{^[R(-),A(-,R(-),W(-)Mm\R(-)}- 
Hence there exists a measurable mapping : C([0, T];R n ) — > R such that 

Z = $(R(-)), Z*=i>{R*{-)). (7.5) 



Proposition 7.2 Let a function <p : C([0,T];R n ) x B([0,T];R n ) x B([0,T];R) -» R 6e suc/i 
£/tcrf E^-(5(-),a(-),r(-)) < +oo. Further, let a function cj> : C([0,T];R n ) xB([0,T];R) K be 
such that E0 _ (i?(-), r(-)) < +oo . T/ien 



E0(i?(-),H(-),r(-)) = EZ*<f>{R*(-),a*(-),n(-)), 
E0(fl(.),r(-)) = E2*0(£*(-),n(-)), 
E^(S(.),r(-))=E^(-),r,(-)). 



(7.6) 
(7.7) 
(7.8) 



Proof. By assumption (O, W(-), a(0)) is independent of To prove (|7.6p it suffices to 

prove 

E^(R(-),a(-),r(-))\@,W(-),a(0)} = E{^0(^(-),a*(-),r*(-))|0^(-)»2(O)} a.s. (7.9) 
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Thus, for the next paragraph, without loss of generality, we shall suppose that (6, W(-), 5(0)) 
is deterministic, since for each value of (G, W(-), 5(0)) we can construct R, R*,a,a*. 

By the linearity of (|2.4p . it follows that Ko(t) defined by (|2.6p is the conditional covariance 
for 5*(i) given (W(-), a(0)). Similarly to the proof of Proposition 12. 1|, it can be shown that (|2.7p 
with m = ensures that E{-Z*|0, W(-), 5(0)} = 1 and Ei?* = 1. We define the probability 
measure P by dP/dP = Z*. (Each value of (0, W(-),5(0)) generates its own P). By Girsanov's 
Theorem, the process 

w(t) = w(i) — / a(s)~ 1 a^(s)ds 
Jo 

is a Wiener process under P. From this we obtain 

dR(t) = A(t, R(-),W(-),a(0))dt + a(t)dw(t), 
dR*(t) = A(t,R*(-), W(-),a(0))dt + a(t)dw(t). 

Then for each value of (O, W(-),5(0)) the processes (R(-),a(-),r(-)) and (i?*(-),5*(-), r*(-)) have 
the same distribution on the probability spaces defined by P and P respectively, and (|7.9p . hence 
(USD follows. 

Further, <\7.7b follows by taking conditional expectation in (|7.6p . Finally, using Proposi- 
tion O and (I7T6D . 

E,^(5(-),r(-)) = E2- 1 0(i?(-),r(-)) = E^(^(-),S(-))- 1 ^(-),r(-)) 

= E^^( J R*(-) J 5*(-))- 1 ^*(-),r*(-)) = E^(i?*(-),r*(-)). □ 

We turn now to Theorem 13.11 Define = F(Z*, A). It follows from (|7.5p that if we define 
by f = then £ = <K#*(0)- 

Proof of Theorem \3.1\ Let us show that EC/ - (£) < oo so that EC/ (£) is well defined. For 
= 1, 2, . . ., we introduce the random events 

n[ k) = {—k < U{t) < 0}, ft (fc) = {—A; < C/(0 < 0}, 

along with their indicator functions, 1**^ and I^ k \ respectively. The number provides the 
unique maximum of the function Z*U(£,*) — A£* over L>, and Xq £ D. By Proposition 17.21 we 
have, for all A; = 1,2,..., 

El( fc )c/(£) - Eli fc) Ae* = Ell fe) fez7(f») - A£*) > Ell fc) (z*C/(X ) - AX ) 
= E#)C/(X ) - XX P(ni k) ) > -\U(X )\ - \XX \ > -oo. 

Furthermore, we have that E|£*| = E*|£| < +oo. Hence EC/~(£) < oo. 
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Now observe that for any n € S we can apply (|7.7p and (|7.8p to U(X 7r (T)) (and use (|7.5p ) 
to obtain 

EU(X*(T)) = E*{ZU(X*(T))} < E*{ZU(X"(T)) - AX^(T)} + XX 

< E m {ZU(£) - A£} + AX = E* ZU(£) = E£/(£). 

Thus (ii) is satisfied. 

Let us show (iii). Since a is non-random, hence w-adapted, then = <fi(w(-)), where 
4>(-) : B([0,T];R") — > R is a measurable functions. By the martingale representation theorem, 

f* = Ef*+ / /(*,<) I [o.t]) T <M*), 

JO 

where f(t,-) : B([0,i];R n ) — > R n is a measurable function such that Jq \f(t,w(-)\[ ^)\ 2 dt < 
+00 a.s. There exists a unique measurable function fo(t,-) : B([0, t];R") — > R n such that 
/(*>™(-)l[o,t]) = /o(i,-R*(-)l[o,t])- Thus, 

£ = E& + f T fo(t,R*(-)\ m ) T dw(t) = E& + [ T f (t,R*(-)\ m ) T a(t)- l dR*(t). 
Jo Jo 

Proposition 17.21 implies that E£* = E*£ = Xq, and 

£=X + F f (t, R(-)\ ^aity'dRit). 
Jo 

Hence the strategy 7r(i) T = B(t)fo(t, i?(-)|[ 0jt ]) T cr(t) _1 replicates B(T)£. It belongs to S; in 
particular, since and R generate the same sigma-algebra and D is convex, then X(t,ir(-)) = 
£ -D, hence bounded below. This completes the proof of Theorem 13.11 □ 

A n(n + l) _ 

Proof of Lemma \4-l\ Let V = R n x R x R 2 . Clearly, V is a n-dimensional linear vector 
space, where n = n + 1 + n(n + l)/2. Let = (yi(t),y 2 (t), y 3 (t)) be a process in V such that 

y(t) = (y(t),y 3 (t)) = (yi(t),y 2 (t)Mt)) = (a(t), hi Z(t),S(t)a(t) T ) . 

The last equality is satisfied by (14. 7p . It can be seen that the equation for y(t) is linear: 

dyi{t) = [^(t)ui(t) + u(t)]cft + E(t) dR(t), 
dy 2 {t) = -\Tr{Q{t)y 3 {t)} dt + yi(t) T Q(t) dfl(t), 

dy 3 (t) = [A(t)y 3 (t) + y 3 (t) T ^) T + v(t)y 2 (t) T + y 2 (^) T + ±{E(t)<r(t)a{t) T E(t) T }]dt 

+E(t) dR(t)y 2 (t) T + y 2 (i) dR(t) T E(t) T . 

Here -A(i), -E(t) are known deterministic functions in R nx ™, R n and R nxri respectively. In 
particular, A(t) = A(t) — b(t)a(t) T Q(t). Thus, the equation for y(t) can be rewritten as 

dm = Km ,t)dt+EtiUm,t)dR t (t), 

y(o) = m, 
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with yo = (mo, 0, mo?rto ), and with some functions f(x,t) : V x [0,T] -> V and gi(x,t) : 
V x [0, T] — * V, i = l,...,n, that are affine in x G V with continuous in t coefficients. In 
particular, df(x,t)/dx and dgi(x,t)/dx depend only on t, and they are uniformly bounded. 
Hence (|7.10p has an unique solution. Therefore, equation (|4.2p has the unique solution y(t). 

Let V(x,t) = E*&(y x,s (T)), where the process y x,s (-) takes values in R n+1 and is such that 
y x ' s (-) = (y x,s (-),y^' s (-)) is the solution of (|7.10p given the initial condition y(s) = x G V. Then 
V(x, t) is the classical solution of the boundary value problem for the corresponding backward 
Kolmogorov's equation 

§(x,t)+C(t)V(x,t) = 0, t G [0, T], 
V(x,T) = $(x 1 ,x 2 ), 

n(n + l) 

where x = (xi,X2,x^) G R n x R x R 2 ; and where C(t) is the second order differential 
operator on functions v : V — > R generated by the Markov process y(t). 

Let y x ' s {-) = (!/?'(■),••• >J/jS+i(0) be the solution of @2J, and let V(x,t) = E*$(y*' t (T)). 
Clearly, 

^'(t) = tf>m,?''®?''(t) T ) = (^(t),...,^ , (t),Ci(*). f' 8 Wf' s (i) T ) » 

if 

X — (x, X3) — (x, X n -)-l , XX ), X — {x\ j • • • , ; — ("El j • • • ; X n , X n _|_l) — (x, X n -|-l ) £ R ^ , 

x 3 = xx T G R«(«+D/2, = (!,*•«(.), . . . , 

In that case, F(x,i) = V"(xi, X2, X2xJ), where x = (x, x n+ i), x G R n . Therefore, V(x,t) is the 
classical solution of problem (|4.3p - (|4.4p . 

Let ?/*(•) denotes the solution of (|4.2p with replaced by -R*(-) = f Q a(t) dw(t). 

Set X*(i) = V(y*(t), i). From (|4.3p and Ito's Lemma, it follows that 

X,(T) = X*(t) + y B,.(s)- l Tr*(s) T dR*(s), 

where vr*(i) T = B*(t)^(y*(t),t)g(y*(t),t). It follows that X*(0) = V(y*(0),0) = 
EV(y*{T),T) = X and 

dX*(t) = B*(t)-\4t) T dR4t), X*(T) = $(y*(T)). (7.12) 

Then X*(t) = ip(t, R*) for some measurable ip, and the result follows if we observe that X(t) = 
ip(t,R) replicates the claim as desired for ir(t) T = B(t)^(y(t),t)g(y(t),t). 
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To continue, we require some a priori estimates. Let C*(0 = B*(t)~ 1 o-(t) T ir*(t). 
We consider the conditional probability space given (B, W(-), a(0)). With respect to the 
conditional probability space, it follows from (|7.12p that 

dX*(t) = Ut) T dw(t), 
X,(T) = <%,(T)). 

By Proposition 2.2 El Karoui et al (1997), the (unique) solution (C*(i), X*(£)) of linear stochastic 
backward equation f)T. 13|) is a process in L2QO, T], L<2(Q, T , P))xC([0, T], L^i^l, F,P)), and there 
exists a constant cq, independent of (<£(•), B, W(-),a(0)) and such that 



su Pte[0jT] E {\X*(t)\ 2 \ B, iy(-), 5(0)} +e|/ T |C*(i)| 



2 dt 



Q,W(-),a(0) 



<c B{<S>(y*(T)) 2 \G,W(-),a(0)} a.s. 

Hence 

sup E|A\(t)| 2 + E / S*(i)- 2 |7r*(i)| 2 di < ciE$(y,(T)) 2 , (7.14) 
te[o,r] Jo 

where ci > is a constant that does not depend on <£(•). Then (|4.5p follows. This completes 
the proof. □ 

Proof of Theorem \4.1\ Clearly, the equation for y(t) is 

da(t) = [A(t)y(t) - b(t)a(t) T Q(t) + a(t)S(t)]dt + ~/{t)Q(t) dR(t), 
dy n+1 (t) = \a{t) T Q{t)a(t)dt -a{t) T Q(t) dR{t). 

As in the proof above, it can be shown that X{t) = V(y(t),t, A) is the solution of some equation 
(|7.12l) . i.e. it is the normalized wealth. Then the proof follows. □ 

Let A/2 be the set of all Gaussian processes a(t) : [0, T] x f2 — > R n which are progres- 
sively measurable with respect to the filtration generated by [a(0),w(t),W(t)] and such that 
E / T \a(t)\ 2 dt < +00. For o(-) G A/" 2 , let 



Z(t, a(-)) = exp 



* 1 /"* 

a{s) T Q{s)dR{s) - - / a(s) T Q(s)a(s)ds 
* Jo 



Proposition 7.3 Xe£ a(-) G A/2, Zei p G (l,+oo), and Ze£ /i G R, /i < or // > 1. Ze£ A"(i) fre 
i/ie covariance matrix of a{t) under P* ; and Zei /c(p) = qT(/j, 2 p — fi), where q = p{p — l) -1 . Let 
n(p)K(t) < o~(t)a(t) T — el n , where e > is a constant. Then E*Z(t, a(-)) M < +00. 
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Proof of Proposition fOl If n G [0,1], then E*Z(i,a(-)) M < +00 ( see Lakner (1998), p.93). 
Therefore, we can assume without loss of generality that fj, < or fj, > 1. Clearly, 



Z(t, a(-)y = exp M /J a(s) T Q(s)dR(s) - % £ a(s) T Q(s)a(s)ds = C(*)Co(*), 



where 



C(t) = exp L [ t a(s) T Q(s)dR(s) - & [ 



Jo 1 Jo 



a(s) Q(s)a(s)ds 



and 




a(s) Q(s)a(s)ds 



By Holder inequality, E*^ < [E*C(T) p ] 1/p [E^oCO 9 ] 17 ". 



Similarly to the proof of Lemma A.l from Lakner (1998), we have that E*£(T) P < +00 
because C(t) p is a positive local martingale with respect to P*, thus by Fatou's lemma it is a 
sup er martingale . 

By Jensen's inequality, 



Kit)- 1 > K(p)[a(t)a(t) T - el n }- 1 = #c(p)0(t)[/ n - eQit)}- 1 

= n(p)Q(t) [I n + Et=i{zQ(t)} k ] > <P)Q{t) + <p)eQ{t) 2 > «(p)Q(t) + M la 



where Mi = M\{e) > is a definitely positive constant matrix. (We can take e > small 
enough to ensure convergency.) Similarly to the proof of Proposition 12.11 it follows from (|7.15|) . 
(pOo]) that E^CoCO 9 < +00 and E,Z(t,a(-)) M < +00. □ 

Proof of Lemma[5J\ If yu G [0,1], then E*Z» < +00 ( see Lakner (1998), p.93). Therefore, 
we can assume without loss of generality that fi < or [i > 1. Note that a(-) G Mi- By 
Proposition 17.31 if (i) is satisfied then E^-Z^ < +00. 

Further, let (ii) be satisfied. Clearly, a(-) G Mi- By Proposition 17.31 again. E*Z(T, a(-))^ < 
+00. By ([73]), J* = E{Z(r,a*(-))|.R*(-),r*(-)}- Hence by Jensen's inequality E*J^ < 
E«Z(T,o(-)y < +00. □ 

Proof of Theorem \6.1\ Let r be a random variable that takes values in [0, T] and such 
that P(r = 0) = 1/2 and P(t G (ii,i 2 ]) = (*a - *l)/(2T) for < ii < t 2 < T. Let r/j G 



E,Co(T) 9 = E.exp \ql&=* d(s) T Q(s)a{s)ds 

= E* exp ^k(p) Jq a(s) T Q(s)d(s)ds <^ E* exp [^K(p)a(s) T Q(s)d(s)] ds. 



(7.15) 



Remind that Q = (co" T ) 1 , and n{p) > 0. Similarly to (|7.ip . we obtain 



(7.16) 
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L<2,(P*, J~, P, R n+2 ) be random vectors such that they have the probability density functions 
Pi(x), i = 0,1. We assume that r, 770, ryi, w, Q, W(-)> a(0) are mutually independent. 
Let 7] = 77oI{ T= o} + ??iI{t>o}i an d let 77*(-) be the solution of the Ito's equation 

dr}*(t) = f(r ] *(t),t)dt + g(r ] *(t),t)dR*(t), t > r, 
t/*(t) = 77. 

Equation (|6.ip is the forward Kolmogorov's equation for the case when time of birth is distributed 
as r, and the vector r]*(t) has the conditional probability density function p(x,t)/2 in the sense 
that P(77*(t) G r, t > t) = 1/2 J r p(x, t)dx for any domain T C R n+1 , where p is the solution of 

fl63p. 

Note that we need random r with the selected probability density on (0, T] to generate the 
free term in parabolic equation (16, ip . 

Assume that $(■) <E C 2 (R n+1 ) and it has finite support. Let V(x,t) = B^(y x ' t (T)), where 
y x ' s (-) is the solution of (|4.2H . Then F(a?,t) is the classical solution of problem (j4.3[) - (j4.4n . Set 
Y*(t) = V(rj*(t), t). From (|4.3p and Ito's Lemma, it follows that 

Y m (T)=Y*(t) + J B*(s)- l Q*(s) T dR*(s), r<t<T, 
where g,(t) T = B,(t)^(^(t),t)g(rj,( y t),t). Hence 

d%{t) = B^ty^^dR^t), n(T) = $(r/*(T)). (7.18) 

To continue, we require some estimates. Let £*(0 = B*(t)~ 1 a(t) T g*(t). 
Consider the conditional probability space given (r, 77, 6, W(-), a(0)). With respect to the 
conditional probability space, it follows from (|7.18p that 

dY*(t) = Ut) T dw(t), 
Y*(T) = HV*(T)). 

By Proposition 2.2 El Karoui et al (1997)) again, the (unique) solution ((*(t),Y*(t)) of stochas- 
tic backward equation UTIM is a process in L 2 ([r, T],L 2 (Q, T, P)) x C([r,T], L 2 (ft, T, P)) 
given (r, 77, 0, W(-), 5(0)), and there exists a constant Co that is independent of 
($(■), T,rj,<d,W(-),a(0)), and such that 



sup te[0jT] EI {t > T} {|n(t)| 2 | t, V, 6, W(-),a(0)} + E|/ T I {t > T} \Ut)\ 2 dt | r, 77, 6, W(-),a(0) 
= sup te[TjT] E {|E(t)| 2 I r,r,, 9, W(-),a(0)} + e| \l(t)\ 2 dt\ r,7?, 9, W(-),5(0)| 

< C E ($MT)) 2 I r, 77, G, ^(0,0(0)} a.s. 
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Hence there exists a constant cq, independent of $(■) and such that 



sup EI {t > T} |y,(t)| 2 + E / l {t > T} BS)~ Z \QS)\ dt < c B$( V *(T)) z . 



(7.20) 



Let $(•) be a general measurable function satisfying the conditions specified in the theorem. 
Then, there exists a sequence where G C 2 (R™ +1 ) are such that they all have 

finite support and 

E|$»(r / ,(r))-$(7 ? ,(T))| 2 = / Rn+2 p(x,r)|$«(x)-$(x)| 2 (ix^0 as i -> oo. (7.21) 
Let (■), £>*(•), and V^(-) be the corresponding processes and functions. Let 



9 itj = sup Ei {t > r} |y (i) (t)-y w (t)| 2 + E / i {t > r} ^(t)^|^(t)-^^(t)| z dt. 

te[0,T] jo 



r 



-2, „(i). 



By (f7T20]) - (f7T2T]) and the linearity of ([713]) . it follows that 

* M <coE|$«(^(T))-$W(^(T))| 2 ^0 as i,j oo. 

We have that G V\£, since they are bounded together with their partial derivatives with 
respect to xi, . . . , x n +i- Remind that < p(x) < p(x,t) for all x,t. Furthermore, we have that 



sup / p{x,t)\V {i \x,t) - V ij \x,t)\ 2 dx 

te[0,T] JR"+! 



fT r 

+ ! dt / p(x,t) 



(a:,i) — 



9a; 



<9x 



g[x,t) 



dx. 



Hence II VW - V^ll* a < tf< -> as z,j — ► oo. Therefore, is a Cauchy sequence in 
Y\}q, and it has the limit V in W^. This V is the desired solution, and (|6.3p is satisfied. This 
completes the proof. □ 

Note that it follows from the proof above that the sequences {y} 1 \-)}°^ 1 and {q*\-)}^i 
are Cauchy sequences in the spaces C([r, T];L 2 (Q, J 7 , P{- | r})) and (["?", £ 2 (^, J 7 , P{- | t})) 
respectively. Hence the corresponding limits Y*(-), £>*(•) exist and belong to these spaces given 

T. 

This paper presents development of some results and ideas that grew up from our collabora- 
tion with Ulrich Haussmann during the author's stay at Pacific Institute for the Mathematical 
Sciences, Vancouver (see, e.g., Dokuchaev and Haussmann (2000)). The author wish to thank 
Prof. U. Haussmann for the support and useful discussion. The author also wishes to thank the 
anonymous referees for their insightful comments which greatly strengthened the paper. 
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